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Abstract 

We consider the amplification of cosmological magnetic fields by gravitational waves as 
it was recently presented in [1], That study confined to infinitely conductive environments, 
arguing that on spatially flat Friedmann backgrounds the gravito-magnetic interaction pro¬ 
ceeds always as if the universe were a perfect conductor. We explain why this claim is not 
correct and then re-examine the Maxwell-Weyl coupling at the limit of ideal magnetohydro¬ 
dynamics. We find that the scales of the main results of [1] were not properly assessed and 
that the incorrect scale assessment has compromised both the physical and the numerical 
results of the paper. This comment aims to clarify these issues on the one hand, while on 
the other it takes a closer look at the gauge-invariance and the nonlinearity of [1], 


1 Introduction 

The interaction between electromagnetic fields and gravitational waves and the possible energy 
transfer between the Weyl and the Maxwell fields have a long research history. A mechanism for 
the amplification of large-scale magnetic fields by gravity waves of similar size soon after inflation 
was recently proposed in [2, 3]. In the poorly conductive environment of early reheating the 
analysis indicated a resonant magnetic amplification proportional to the square of the field’s scale 
(see [3] for details). This meant that Weyl-curvature distortions could provide a very efficient 
early-universe dynamo of superhorizon-sized magnetic fields. For example, fields with a current 
comoving scale of approximately 10 kpc and a strength of ~ 10~ 34 G, like those produced in [4], 
could be amplified by many orders of magnitude by the end of reheating. 

The same gravito-magnetic interaction has been applied to infinitely conductive cosmologies 
in [1], Central to that study is the claim that on spatially flat Friedmann-Robertson-Walker 
(FRW) backgrounds the Maxwell-Weyl coupling proceeds always as if the universe were a perfect 
conductor. We argue that this is not the case and explain why the above mentioned paper arrived 
at the opposite result. In addition, we find serious problems in the scale assessments of [1]. 
These have prevented the authors from recognising the main effect of the interaction, namely 
the large-scale superadiabatic amplification of the R-field, and compromised their results. With 
the present comment we draw attention to these issues by correcting the mathematics, where 
necessary, and by clarifying the physics. We also re-examine and question the gauge-invariance 
and the nonlinearity of the formalism proposed in [1], 
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2 On the electric curl and the role of conductivity 


Section III.B of [1] argues that the Maxwell-Weyl coupling proceeds unaffected by the conduc¬ 
tivity of the medium and as if the universe were a perfect conductor. This follows the claim that 
the gravitationally induced, second order, electric field ( E a ) is curl-free if the zero order FRW 
model has flat spatial sections. This is not the case, however, because 

(curl. E a )' = -0curlE a + lZ ab B b - curl J a - D 2 B a , (1) 

at second order. Here 0 is the volume expansion, lZ ab is the first-order 3-Ricci tensor, B a is 
the first order magnetic field, J a is the spatial current and D 2 = D a D a is the 3-dimensional 
Laplacian operator. Therefore, even when curlE a initially vanishes, there are sources in the 
right-hand side of (1) that will generally lead to a nonzero electric curl. Switching the latter off 
is not a consistent constraint. Note that the current term depends on the conductivity of the 
medium, while the 3-curvature distortions and the magnetic field fluctuations are caused by the 
presence of gravity-wave perturbations. These set the field lines into motion causing magnetic 
fluctuations that produce an electric component. The latter has curlE a / 0 because of (1). This 
means that the gravito-magnetic interaction does not always proceed as if the conductivity of 
the universe were infinite. In [1], the authors arrived at the opposite conclusion because they 
considered the second time-derivative of curl£' a instead of the first, namely (see Eq. (29) in [1]) 

(curlE a )" + |0(curl£’ a )' - D 2 curl.E a + [|0 2 + - 9 p) + |A] curl. E a = curlA' a . (2) 

In the above, which is said to hold at the second perturbative level, the pair /i and p represents 
the density and the pressure of the matter, A is the cosmological constant and K a is a gravito- 
magnetic source term. Arguing that curliv a vanishes when the FRW background is spatially 
flat (see Eq. (30) in [1]), the paper claims that the electric field will remain curl-free if it was 
so initially. However, this a priori sets (curlE a )' = 0 in Eq. (2), although the latter is not 
necessarily the case because of expression (1). 

By switching the electric curl off the authors have inadvertently confined their study to 
idealised perfectly conductive universes, bypassing all plasma effects and the implications of finite 
conductivity. Because of that [1] cannot follow the gravito-magnetic interaction in the poorly 
conductive stages of early reheating and therefore the comparison with [2, 3] was inappropriate. 


3 On the scale assessment and the nature of the amplification 


Restricting to perfectly conducting cosmological environments we have (see § II.C, IV.B in 

B a + §©-B a = (J ab B b = la , 


[ 1 ]) 

( 3 ) 


where B a is the total magnetic field, B a is the original and a a b is the transverse shear component. 
Then, the gravitationally induced B-fields during the radiation era is 


Br = Bo 









( 4 ) 
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while for dust and late reheating we have 


Bd/rh — Bo 




ClQ \ 3 / 2 

a ) 
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(see Eqs. (46), (47) in [1]). Above a is the scale factor, H is the Hubble parameter and the zero 
suffix marks the onset of the gravito-magnetic interaction. These expressions were treated as 
infinite wavelength solutions (see § IV.B.l in [1]), when in reality they cover all finite superhorizon 
scales. 1 This incorrect scale assessment has led to the sequence of problems that we will discuss 
below. Here we note that the domains of (4) and (5) are essential for the additional reason that 
both results show a deviation from the standard a~ 2 magnetic decay-law and a superadiabatic 
amplification of the field (i.e. B = constant for radiation and B tx a -1 for dust/reheating). 

Returning to [1] we find (see Eq. (51) there) that the final solution for the gravitationally 
induced magnetic field, during both the radiation and the dust eras, is 



1 
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since p 7 oc a” 4 is the radiation density. 2 It is also argued (see § IV.B.2 after Eq. (51) in [1]) 
that the above spans all finite wavelengths irrespective of their size relative to the horizon. This 
cannot be correct because (8) means B oc a~ 2 on all scales, while (4) and (5) give B = constant 
and B oc a~ 1 on super-Hubble lengths. Consider, for example, a large-scale B-mode propagating 
through the radiation era and ignore any effects from reheating. Then, Eq. (4) gives 


B ~ Bo 


1 + 
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Well outside the horizon Xh/X^ -C 1, which means that using (8) instead of (9) on those 
scales will grossly overestimate the amplification of the field. The two expressions agree only at 
horizon crossing when Xg = Xh- Therefore, at best, Eq. (8) works for scales that have crossed 
the Hubble radius. However, the magnetic growth takes place earlier outside the horizon and 
results from the superadiabatic amplification of the field seen in (4) and (5). In [1] the reader 

1 Assuming dust and setting r = t/to the gravito-magnetic source term in the right-hand side of (4) reads 

I W (t) = [lb J 5/2 (2fr 1/3 /a 0j ffo) + B 2 Y 5/2 (2£r 1/3 /aoB 0 )] r" 5/2 , ( 6 ) 

where l the wavenumber of the gravitational wave (see Eq. (48) in [1]). Finite superhorizon scales have 
It 1 / 3 /do-ffo -C 1 and the series expansion of the Bessel functions (with the initial conditions of [1]) give 

/(t) = 2aoBor~ 5 ^ 3 — aoBor^ 10 ^ 3 . (7) 

Inserted into Eq. (4) the above leads to solution (5), which therefore applies to all finite super-Hubble lengths 
and not only to infinite wavelengths. Similarly one can show that (5) also covers all finite superhorizon lengths. 

2 Expression (8) was obtained after replacing A gw with A^ in Eqs. (49) and (50) of [1]. However, this substitu¬ 
tion was arbitrary because no relation between the two scales was previously established and the original magnetic 
field was given a zero wavenumber (see § IV.B in [1]). The latter means that A^ is ill defined (i.e. A^ > 00 ) and 

misleadingly suggests an infinitely strong induced B-field in (8). Note that assigning a nonzero wavenumber to 
B provides a useful relation between As, A gw and A g (see Eq. (10) in [3]). 
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is unaware of the nature of the magnetic growth and where it occurs. The root of the problem 
was assigning (4) and (5) to infinite wavelengths. As a result, these solutions were sidestepped 
(see in particular § V after Eq. (60) in [1]), the suparadiabatic nature of the amplification was 
not recognised, the domain of (8) was incorrectly assessed and ultimately the physics of the 
gravito-magnetic interaction was misrepresented. 

An additional complication is that Eq. (8) alone cannot provide the total amplification of a 
mode that went through different epochs (e.g. reheating and radiation), because each domain 
has different growth rates (compare (4) and (5)). Consider an inflationary H-mode that crosses 
the horizon in the radiation era (like that in § V of [1]). Assuming for simplicity that the growth 
is always strong, the residual field according to (8) is B = Bo(cr/H)o(X^/ A#)g(ao/a) 2 . This 
result was used in [1] but it does not incorporate the reheating effects, for which we also need 
(5). Then, 


B = Bn 


-) 

H) 


H ) RH 


(&)!( 


aoV 

a ) 


( 10 ) 


where RH marks the end of reheating. The difference in the above is made by extra factor 
(a/H)RH. The latter is typically very small and this accounts for the fact that the magnetic 
growth rate during (late) reheating is considerably slower than that of the radiation era. There¬ 
fore, not appreciating the role of (4) and (5) has also compromised the numerical results of [1], 


4 On the gauge-invariance and the nonlinearity 

Studies of cosmological perturbations are known to suffer from the gauge problem. The aim 
of [1] is to provide a nonlinear treatment of the gravito-magnetic interaction free from gauge 
ambiguities. This meant integrating the magnetic induction equation (i.e. Eq. (3) above) with 
respect to B a (see § IV in [1]). However, B a is treated as the perturbation of B a (see § II.C 
in [1]). The latter has nonzero linear value and this makes B a a gauge-dependent vector at 
second order by known theorems [6]. In an attempt to circumvent the problem the auxiliary 
variable (3 a , with (3 a = B a + 2QB a /3 was temporarily introduced (see § II.C in [1]). This has 
zero linear value and is therefore gauge-invariant at second order. Nevertheless, (3 a has been of 
little practical use because the authors still had to solve for the gauge-dependent vector B a to 
extract any meaningful information about the evolution of the field (see Eqs. (4), (5) and (8) 
here). 3 This fact makes the analysis and the results of [1] gauge dependent. 

The nonlinearity of [1] is built on a set of four spacetimes (see § II there). However, the new 
setting has not improved our understanding of the interaction because still only the gravito- 
magnetic effects of [2, 3] are accounted for and all other nonlinearities (including the magnetic 
effects on the shear) are excluded. There is no new information in the nonlinear equations 
of [1], relative to what is already encoded in the linear formulae of [2, 3]. 4 This should not have 

3 If the authors had found a way of evaluating the magnetic growth by means of f3 a exclusively, their results 
would have been gauge-invariant. This was not possible however. Moreover, at the ideal magnetohydrodynamics 
(MHD) limit f3 a is all but redundant because f3 a = &abB b (see § II.C, IV.B in [1] or Eq. (3) here). 

Expressions (4) and (10), the latter with A = 0, of [1] are identical to Eqs. (2b) and (4) in [3]. Also, at the 
MHD limit, the magnetic wave equation of [3] (see Eq. (3) there) and expression (12) in [1] reduce to Eq. (3) here. 
Moreover, the constraints used in [1] (see § II.B.2 there) are the standard linear ones. These constraints ensure 
that the traceless tensors remain transverse at all times, a highly nontrivial issue for any truly nonlinear study. 
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happened and the reason it does is the selective setting of [1], which restricts the study and 
prevents it from going beyond the linear level. 

5 Discussion 

The Maxwell-Weyl coupling and the possible energy transfer between the two fields has a long 
research history. Recently this interaction was proposed as a very efficient (resonant) ampli¬ 
fication mechanism of large-scale magnetic fields during the poorly conductive stages of early 
reheating [2, 3]. The same mechanism was studied at the ideal MHD limit in [1], claiming that 
the gravito-magnetic interaction proceeds always as if the universe were a perfect conductor 
when the background spatial geometry is Euclidean. We explained here why the aforementioned 
claim is not correct and how it restricts the generality of [1], We also demonstrated that the 
aforementioned paper did not properly monitor the gravito-magnetic interaction. In particular, 
solutions which in [1] were assigned to infinite wavelengths only, were found to hold on all finite 
super-Hubble scales. Also solutions that work only inside the horizon were applied to all finite 
scales. As a result, the physical interpretation and the numerical results of [1] were seriously 
compromised. With this comment we have attempted to clarify these issues by correcting, where 
necessary, the mathematics of the analysis and by explaining the physics of the interaction. In 
the process we also considered and questioned the gauge-invariance and the nonlinearity of the 
formalism proposed in [1], 
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